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Search for the proton decay p K^v in the large hquid scintillator low energy 

neutrino astronomy detector LENA 

T. Marrodan Undagoitia/'E F. von Feilitzsch,i M. Goger-Neff/ C. 
Grieb,^ K. A. Hochmuth/ L. Oberauer/ W. Potzel,^ and M. Wurm^ 

Physik- Department E15, Technische Universitdt Miinchen, 
James-Franck-Str., 85748 Garching, Germany 
"^Virginia Polytechnic Institute and State University, 
Physics Department, Robeson Hall, Blacksburg, 
VA 24O6I-O435, USA 

The LENA (Low Energy Neutrino Astronomy) detector is proposed to be a large- volume liquid- 
scintillator device which will be highly suitable for the investigation of a variety of topics in astro- 
physics, geophysics and particle physics. In this paper, the potential of such a detector concerning 
the search for proton decay in the SUSY favored decay channel p — » K^V is investigated. Based on 
Geant4, Monte Carlo simulations of the proton decay in the LENA detector as well as of the back- 
ground radiation in the detection energy windows have been developed. From these simulations an 
efficiency of 65% for the detection of a possible proton decay has been determined. Within ten years 
of measuring time a lower limit for the proton lifetime, concerning the decay channel investigated, 
of r > 4 • 10'^'' y could be reached. 
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I. DETECTOR DESCRIPTION 

The LENA (Low Energy Neutrino Astronomy) detec- 
tor is planned to have a cylindrical shape with about 
100 m length and 30 m diameter. An inside part of 13 m 
radius will contain approximately 50 kt of liquid scin- 
tillator while the outside part will be filled with water 
to act as a muon veto. A fiducial volume for proton 
decay will be defined having a radius of 12 m. Cover- 
ing about 30% of the surface, 12 000 photomultipliers of 
50 cm diameter each will collect the light produced by 
the scintillator. PXE (phenyl-o-xylylethane) is foreseen 
as scintillator solvent because of its high light yield and 
its save handling procedures. The optical properties of a 
liquid scintillator based on PXE have been investigated 
in the Counting Test Facihty (CTF) for BOREXINO at 
the Gran Sasso underground laboratory ^2]. A yield of 
372 ± 8 photoelectrons per MeV (pe/MeV) have been 
measured in this experiment with an optical coverage of 
20%. The attenuation length of ^ 3 m (at 430 nm) was 
substantially increased to ~ 12 m purging the liquid in 
a weak acidic alumina column Q . With these values an 
expected photoelectron yield of ~ 120 pe/MeV can be 
estimated for events in the center of the LENA detector. 
Currently the optical properties of mixtures of PXE and 
derivatives of mineral oils are under investigation |^. 

In addition to astrophysical processes such as super- 
novae neutrinos, relic supernovae neutrinos and solar 
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neutrinos also geophysical aspects and research in parti- 
cle physics, e.g. long baseline experiments and the search 
for the proton decay will play a major role. Possible loca- 
tions for this detector, currently under discussion, are an 
underground mine in the center of Finland (Pyhasalmi, 
CUPP: Center of Underground Physics in Pyhasalmi) or 
under water, close to the coast at Pylos in Greece. Both 
sites have a shielding of ~ 4 000 meters water equivalent 
(m.w.e.) and are far away from nuclear power plants, 
which cause the main contribution to the Vf, background 
in the measurement of relic supernovae neutrinos. 

In this letter the sensitivity of LENA for the search of 
the proton decay p — > K^V is investigated. The potential 
of a large liquid scintillator detector for proton decay 
detection has already been proposed by R. Svoboda . 



II. THEORETICAL PREDICTIONS 

In the Grand Unified Theories (GUTs), the decay of 
the proton is predicted, leading to a nonconservation of 
the baryon (B) and lepton (L) numbers. A measure- 
ment of such an event could be a probe of the validity 
of those theories and further evidence for the existence 
of physics beyond the Standard Model. In the minimal 
GUT SU(5) [g, the nucleon decays via exchange of gauge 
bosons with GUT scale masses (~ 10"'^^ GeV). The dom- 
inant proton decay mode predicted is p — > e^ir^ with a 
lifetime of t 10"^^ y. The Super-Kamiokande Collabo- 
ration has set the current lower limit for this channel to 
T > 5.4-1033 y at 90% confidence level (C.L.) This re- 
sult together with the fact that, in this theory, the three 
running coupling constants do not meet in a single point 
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give reasons to look also at other theories. 

In Supersymmetric GUTs (SUSY), the nucleon de- 
cay is affected by operators of dimension 4, 5 and 6. 
They involve B- and L-violating processes but conserve 
B-L. For the proton decay mode p — > e~''7r'^ a lifetime 
of r ^ 5 • 10'^^*^ y has been calculated. This rather 
long lifetime results from the large value of the unifica- 
tion mass, Mq, given by those theories. For the decay 
channel p — > K^V {p = "Yl, {i=e t}^) ^ range for the 
lifetime r ~ (0.33 — 3) ■ 10'^^ y has been predicted. 
Also previous studies on SUSY theory have found the 
same decay channel p — > A'+F to be the dominant one 
with a lifetime range of r ^ (10^° - 10^^ y) The kaon 
produced in the dominant SUSY proton decay channel is 
'invisible' in water Cherenkov detectors because its mo- 
mentum is below the threshold for producing Cherenkov 
light {Tth{¥.+) = 253 MeV). For this reason, the Super- 
Kamiokande experiment could only perform an indirect 
measurement with efficiencies of 4.4% and 6.5% corre- 
sponding to two different methods Q. The lower limit 
that they reach at 90% C.L. is r > 2.3 ■ 10^3 y jl3|. 

Supergravity theories provide also predictions for the 
possible decay channels of the proton. The dominant 
mode here is p — > tt"*"!/ with a branching ratio of 65.7%. 
However, again the proton decay mode p — > K'^V is pre- 
dicted with a branching ratio of 33.5% 

In a scintillation detector like LENA, the K"*" particles 
from proton decay events can be observed directly. These 
signals together with the delayed signals from the decay 
products provide a very distinct signature for this decay 
mode. 



III. DETECTION MECHANISM 

Within the detector volume of LENA, about 
1.45- 10'^'* protons, both from C- and H-nuclei, are can- 
didates for the decay. This number has been calculated 
taking into account the chemical composition of PXE 
(CigHis), its density (0.985 g/cir?) and its molecular 
weight (210.3 a.m.u.). 

In the case of protons from H-nuclei (0.23 ■ 10'^^ pro- 
tons in the fiducial volume of LENA), the decay problem 
is simplified because those protons decay at rest. There- 
fore, the proton decaying in the channel p K^V can 
be considered as a two-body decay problem where 
and V always receive the same energy. The energy cor- 
responding to the mass of the proton, rup = 938.3 MeV 
is thereby given to the decay products. Using relativistic 
kinematics, it can be calculated that the particles receive 
fixed kinetic energies, the neutrino 339 MeV and the kaon 
105 MeV. 

In the LENA detector, the kaon will cause a prompt 
mono-energetic signal while the neutrino escapes without 
producing any detectable signal. After t{K^) = 12.8 ns, 
the kaon decays via — > jJL^v^ with a probability of 
63.43 %. In 90% of the times, the kaon decays at rest 
and from this kaon decay, a second mono-energetic sig- 



nal arises, corresponding to the /i"*" with 152 MeV kinetic 
energy. Later, after T(/i+) = 2.2 /is, also the muon will 
decay: p,^ e^VeV^^, the e+ will produce a third long- 
delay signal. Since the last decay is a three body decay, 
the e"*" does not have a fixed energy but a spectrum. 
With a smaller probability, 21.13 %, the kaon decays, al- 
ternatively, via — > tt+tt". The 7r+ deposits its kinetic 
energy (108 MeV) in the detector. The tt" decays al- 
most immediately, within t{tt^) — 8.4 ■ 10^^ ns, into two 
gammas with a total energy being equal to the sum of 
the kinetic energy of the 7r° (111 MeV) and its rest mass 
(tottO = 135 MeV). These gammas can be seen in the de- 
tector as electromagnetic showers. Thus again there are 
two mono-energetic signals, firstly the kaon and secondly 
the one corresponding to the sum of the kinetic energy 
of the TT^ and the energy of the gammas. Again, there 
will be a third long-delayed signal from the succeeding 
decays (tt^ — > ^'^v^ and /j,"*" — > e'^VeV^). Only these two 
most likely decay channels of the kaon are considered in 
the analysis. 

For the protons of the C-nuclei (1.22 • lO"^** protons 
in the fiducial volume of the LENA detector), one has 
to consider further nuclear effects on the proton decay. 
First of all, since the protons are bound to the nucleus 
their effective mass will be reduced by the binding energy. 
The proton energy available for distribution between the 
decay products, will be iW^ — Mp — Ei,, where is 
the modified proton mass, Mp — 938.3 MeV is the rest 
mass of the proton and Ef, is the nuclear binding energy. 
The value of the binding energy Et can be derived from 
gaussian probability functions which for C-nuclei ^3 are 
centered at 37 MeV and 16 MeV for protons in s-state 
and p-states, respectively. Secondly, in case of a proton 
decay in a nucleus, e.g., carbon, decay kinematics are 
different from the free proton due to the Fermi motion 
of the proton. These Fermi momenta have been mea- 
sured by electron scattering on ^^C [l^ . The maximum 
momentum is about 250 MeV/c. 

To determine the kinetic energy change of the decay 
particles of the proton, various decay processes have been 
calculated First of all, the mass of the proton was 
modified due to the binding energies for s- and p-states 
according to the values mentioned above. For both, two 
extreme cases were considered: the neutrino, after the 
decay, moving in the direction of the original proton and 
the neutrino moving in the opposite direction, the kaon, 
of course, always moving in the opposite direction of 
the neutrino. The momentum of the proton was cho- 
sen to be maximal, corresponding to a Fermi momentum 
of ~250 MeV/c. The limiting values for the range of the 
kinetic energy of the kaon are: (25.1 — 198.8) MeV for 
the s-state and (30.0 - 207.2) MeV for the p-state. 



IV. SIMULATION 

In order to quantitatively estimate the potential of 
the LENA detector for measuring the proton lifetime, a 
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Monte Carlo simulation for the decay channel p K'^v 
has been performed. For this pur pose, the Geant4 sim- 
ulation toolkit has been used The event generator 
was programmed to produce K"*" particles of 105 MeV in 
the center of the detector volume and with random di- 
rections. The simulation was also performed for proton 
decays which do not take place in the center to investi- 
gate the dependence on the event position. 

Not only all default Geant4 physics lists were included 
but also optical processes as scintillation, Cherenkov light 
production, Rayleigh scattering and light absorption. For 
the last two processes scattering and absorption lengths 
of As = 60 m and Aq — 12 m, respectively, were assumed. 
With these numbers an effective attenuation length of 
A = 10 m results from: 



This is a quite conservative approach compared to the 
measured value of A = 12 m Q and results in a light yield 
of ~ 110 pe/MeV for an event in the center of the detec- 
tor. The dependence of the detection efficiency and the 
photoelectron yield on scattering and absorption lengths 
will be discussed later (see Chapter 6, Table HTjl . 

To determine the detector efficiency, a total of 20 000 
p K^V Monte Carlo events were generated. Within 
this simulation, it was possible to verify that the kaons 
decay according to the branching ratios for the different 
channels stated in J5,J . The amount of light emitted by 
a scintillating material is not strictly proportional to the 
energy deposited by the ionizing particle. In reality, it is 
a complex function of energy, the type of particle and the 
specific ionization . To take into account the so called 
quenching effects, the semi-empirical Birk's formula [l7| | 
has been introduced into the code: 



V. BACKGROUND REJECTION 

A. Muon Production 

The main background source in the energy range of 
the proton decay are muon neutrinos produced by 
cosmic ray interactions in the atmosphere. Atmospheric 
can interact with the scintillator producing muons in 
the energy range where the search for the proton decay 
is performed. The rate of these events in the relevant 
energy range can be derived from the Super-Kamiokande 
measurements 0- From this data we find an event rate 
of: 

F - 4.8 • 10"2(MeV"^kt" V"^) (3) 

where the LENA energy window as well as the volume 
have to be introduced in MeV and kt, respectively. 

In order to distinguish the real proton decay signals 
from muon background events, a pulse shape analysis 
can be applied. From the simulation it is known that the 
kaon deposits its energy in 1.2 ns. Then, as the kaon de- 
cays within a relatively short time (t(K~'') = 12.8 ns), the 
signature from such an event will consist of two shortly 
delayed signals. Examples of a proton decay signal where 
the separation is clearly visible, where it is almost hidden 
in the risetime and of a muon background signal can be 
seen in figure The better the two signals can be sepa- 
rated, the better the background suppression will be. In 
the cases where the kaon decays at late times (t > 10 ns), 
the two peaks can be easily separated; in all other cases, 
only a shoulder within the risetime can be seen. 

The on-set time of these examples is 60 ns as the events 
are generated at the center of the detector and the pho- 
tons need this time to reach the photomultipliers. 
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The parameters of this formula are A: absolute scin- 
tillation efficiency and kB: parameter relating the den- 
sity of ionization centers to ^ . Both are specific values 
for every scintillator and are obtained from experimental 
data. For the simulation, A = photons/MeV and 
kB = 1.5 • 10~^ cm/MeV have been taken from studies 
performed by the BOREXINO Collaboration Ac- 
cording to equation |21 the program calculates in every 
step the deposited energy (dE) within the distance {dx) 
travelled by the particle and derives the light output 

fdL\ 

Also the photomultipliers have been introduced in the 
simulation assuming an efficiency of 17 % and a time 
jitter of (T = 1 ns originating from the transit time of 
electrons through the photomultipliers and from the sta- 
tistical noise. These features are known because this type 
of photomultiplier has already been used in other exper- 
iments |l9f . 



1. Time Cut 

A first step towards a pulse shape analysis has been 
taken in order to quantify the signal and background dis- 
crimination. Not only 20 000 Monte Carlo proton decay 
events were generated but also 20 000 muon background 
events for each kaon decay channel analysed. For these 
simulations, the particles were generated in the center of 
the detector and the light produced was followed through 
the scintillator and the photomultipliers. This results in 
the number of photoelectrons detected as a function of 
time. From the resulting signal the maximum number 
of photoelectrons is taken and the 15 % and the 85 % 
values of the maximum signal height are calculated. The 
difference of the two corresponding time values is plotted 
in figure 121 

As an example, a cut at 7 ns is shown in this figure. 
Applying such a cut, a background reduction by a factor 
of ^ 2 • 10* is achieved. The background plot shown 
applies only to the channel p — > K^V and —f ^.^v^. 
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A similar simulation has been performed for the channel 
p —> if +17 and — > tt'^tt'^ resulting in an equivalent 
result where the same rejection factor is obtained. An 
efficiency of St = 0.65 results after this time cut because 
only 35% of the signals are lost. 
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FIG. 1: Proton decay signal where the kaon decays after 18 ns 
(top) and after 5 ns (middle). The bottom panel shows a 
background signal. 



2. Energy Cut 

In order to further improve the signal to background 
separation a cut in the energy spectrum can be applied. 
Figure 13 shows the distribution of the number of pho- 
toelectrons coming from free proton decay events at the 
center of the detector as seen in the photomultipliers. 

Two peaks can clearly be identified, corresponding to 
the two main decay channels of the kaon. They are re- 
lated to the energies of 257 and 459 MeV, which is the 
sum of the energies deposited by the kaon and its decay 
products (see section Detection Mechanism). For each 
peak a 300 MeV window has been set resulting in an 
efficiency of e^; = 0.995 for this energy cut. The two 
energy windows overlap resulting in a final window of 
500 MeV. The signals outside these energy windows cor- 
respond mainly to kaon decay channels not considered in 
the analysis. There is also a small loss due to inelastic 
hadronic interactions of the decay particles before they 
are stopped. 
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FIG. 2: Distribution of time intervals in the risetimes for 
proton decay (top) and muon background signals (bottom). A 
cut at 7 ns yields a reduction factor of ~ 2 ■ 10^ on background 
events whereas the efficiency for observing a proton decay 
p K+V is 65%. 
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B. Hadron Production 

Atmospheric neutrinos can interact with the detector 
producing also hadrons. The most probable of these re- 
actions is the single pion production [23| [l^l : 



(4) 



In such a reaction, the first signal comes from the sum 
of the energies deposited by the /i^ and the 7r+. Later, af- 
ter = 26 ns, the tt"*" decays into /i^ and one neutrino 
producing a second short-delay signal. However, the tt^" 
has such a small mass (139.6 MeV) that its decay muon 
receives only about 20 MeV kinetic energy. The signal 
produced by the /x"*" is small. Therefore, those events do 
not have the same signature as the proton decay signal. 
Moreover, since the muon signal is so small it is usually 
hidden behind the above mentioned first signal. 

Neutrinos of higher energies can interact producing 
strange particles, thus also kaon production has to be 
considered. According to the MINERvA proposal [2^ . 
the reactions of this type with notable probability are: 
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FIG. 3: Distribution of the number of photoelectrons pro- 
duced in the photomultipliers coming from free proton de- 
cay events at the center of the detector. The cuts performed 
for the two main decay channels of the kaon (at ~ 257 and 
~ 459 MeV) are shown. It can be recognized that the en- 
ergy windows taken overlap resulting in a total window of 
55000 photoelectrons ~ 500 MeV. 



In equations|Bland[71 to conserve the strangeness num- 
ber the kaon is produced together with a A° baryon. This 
A" particle decays after r^o — 0.26 ns mainly via two 
channels A" ^ p + tt" (63.9%) or A" ^ n + 77° (35.8%). 
Because of this short lifetime, the p, tt^ or n, 7r° cause a 
prompt signal together with the the and the vr" 
in the case of equation[71 The prompt signal in these two 
reactions is bigger than the signal of the /i+ from the kaon 
decay (or 7r+, tt" if the kaon decays within the second- 
probable channel). For this reason, the signature of the 
background signal can be distinguished from proton de- 
cay events. Figure 01 shows an example of equation |S1 
First a big peak can be seen corresponding to the signals 
of the /.t", and A". After that, two other peaks ap- 
pear, one is due to the ^.^ that originates from the decay 
of the if + and the other a fjL~ from hP — > p -|- 7r~ and 
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FIG. 4: Background signal for proton decay events. The first 
peak corresponds to the signals of the fj,~ , and A". At 
later times, two other peaks appear, one is caused by the jj.'^ 
that comes from the decay of the K'^ and the other a fj.~ from 
A° ^ p + Tv^ and tt" /i^F,!. 

Reaction can be responsible for a potential back- 
ground for proton decay in the LENA detector. For 
this case, only neutrinos with energies between 650 and 
900 MeV can produce a signal with a signature similar 
to that of the proton decay. This neutrino energy win- 
dow has been determined considering the energy win- 
dow of proton decay events relevant for the analysis 
(150 — 650 MeV) and the energies deposited by the parti- 
cles in the possible background events. The atmospheric 
neutrino flux has a maximum around 100 MeV and then 
decreases exponentially to higher energies. Assuming an 
energy dependence of 0^ ^ E~'^-'^ for the neutrino flux 
and a linear dependence for the cross-section, 10% of the 
total atmospheric neutrino flux occurs within the energy 
window considered in our analysis. Within this window 
a rate of 0.8 events per year {y~^) caused by such neu- 
trino reactions is predicted. To distinguish these back- 
ground reactions from proton decay events the number 
of delayed electrons produced will be taken into account. 
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For proton decay in the channel considered, always one 
and only one electron is produced by the decay of the 
kaon. In the background reaction always two electrons 
are present, one from the /i~ decay and another from 
the decay chain. This fact considerably reduces the 
background contribution. Only in those cases where the 
energy of one of the background electrons is very small 
(E < 0.5 MeV) or when the signal of these electrons is 
hidden under that of the parent particle, the process ex- 
hibits only one electron. Events with electrons being 
hidden within the previous signal happen in 4% of the 
cases. From 200 muons simulated none of them had a 
decay electron with an energy smaller than 0.5 MeV. To 
estimate the number of background events of this type, 
the rate predicted has to be multiplied by this 0.04 and 
then, by two because electrons both from the fi^ and the 
K'^ decay chain can be hidden within the previous sig- 
nal. Thus, the number of events of the type of reaction^ 
expected in the LENA detector is 0.064 y~^. 

Other reactions where hadrons are produced are pos- 
sible but all of them have lower rates and can also be 
distinguished from a proton decay signal by applying the 
arguments used for reactions El and [7| 



C. Neutron Production by Cosmic Muons 

Also very fast neutrons, generated by spallation pro- 
cesses of cosmic muons outside the detector, may reach 
the scintillator causing background. These neutrons are 
slowed down through scattering processes and the recoil 
protons may produce a background signal in the detec- 
tor. However, only neutrons produced by muons passing 
close to the detector but not through the veto, have to 
be taken into account. In addition, these signals can 
be tagged as the neutrons are captured by hydrogen in 
the scintillator after ^ 200 /xs, producing an observable 
2.2 MeV 7-ray 



VI. PROTON DECAY SENSITIVITY 

Once the time cut and the energy cut have been ap- 
plied, one can calculate the final efficiency e — Ee ■ £t of 
the proton decay detection in the channel considered. In 
tablelHthe efficiencies for events at different points of the 
detector are presented as well as the background rates 
for the same positions. As the detector is symmetrical 
along the cylinder axis the points considered are in the 
middle of the cylinder and at different radial distances. 
No dependence of the efhciency on the event position 
was found within the statistical error. This statistical er- 
ror mainly comes from the time cut as the one from the 
energy cut is negligible. The background rates presented 
have been obtained by multiplying the muon background 
rate in the LENA detector (1190.4 y^^) derived from the 
Super-Kamiokande experiment |9| with the background 
suppression from the time cut (5 ■ 10~^). The contribu- 



TABLE I: For 1000 events produced at each radial position 
R in the detector, the efficiency in the time-cut et, the ef- 
ficiency in the energy-cut ee, the total efficiency e including 
the statistical error, and the background rate B are given. 



R (m) 


Et 


EE 


E (±0.04) 




Center 


0.649 


0.995 


0.65 


0.11 


3.0 


0.675 


0.994 


0.67 


0.11 


6.0 


0.650 


0.994 


0.65 


0.11 


9.0 


0.679 


0.995 


0.68 


0.11 


11.5 


0.666 


0.995 


0.66 


0.11 



TABLE IL Efficiencies e for the detection of proton decay 
and photoelectron yield as function of attenuation, absorp- 
tion and scattering lengths. The last column shows the time 
(see figure 121 at which the cut can be performed such that all 
background signals are rejected. Upper part: Different atten- 
uation lengths for Aa = A^. Bottom part: Attenuation length 
A=10 m for different combinations of Aa and As. 



A (m) 


Aa (m) 


As (m) 


E 


Y (pe/MeV) 


Cut (ns) 


5 


10 


10 


0.56 


58 


10 


7 


14 


14 


0.65 


116 


8 


9 


18 


18 


0.67 


161 


7 


10 


12 


60 


0.65 


110 


7 


10 


15 


30 


0.69 


145 


7 


10 


20 


20 


0.66 


180 


7 


10 


30 


15 


0.63 


230 


8 


10 


60 


12 


0.62 


303 


9 



tion of the kaon production background (0.064 y ^) has 
been added to the previous value. 

In table ^1 the efficiencies for proton decay detection 
and the photoelectron yield are given as function of differ- 
ent attenuation, absorption and scattering lengths. For 
the results presented, 1000 proton decay and 1000 back- 
ground events were simulated for each combination of Aa , 
Xs and A. All events were programmed to take place at 
the center of the detector. The last column of the table 
shows the time (see figure at which the cut can be 
performed such that all background signals are rejected. 
In the upper part of the table, three attenuation lengths 
A are presented keeping the contributions of absorption 
length and scattering length equal, Xa ~ Xg. As the at- 
tenuation length increases the photoelectron yield Y in- 
creases and the cut in the time intervals can be performed 
earlier. In the bottom part, the attenuation length has 
always the same value A = 10 m and different values 
for the absorption and scattering lengths are given sat- 
isfying equation (1). The photoelectron yield increases 
as the absorption length increases because more optical 
photons reach the photomultipliers and therefore more 
photoelectrons are produced. As the scattering length 
increases, the cut in the time intervals can be performed 
earlier because the time information is better. 

Table Hi) shows that in all cases investigated, for A > 
7 m the efficiencies e > 0.62. This is also true for 
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Xa = 12 m and = 60 m, which we used for de- 
riving the resuhs presented earher in figures ^ and |21 
It demonstrates that the efficiency reachable is about 
tenfold enhanced compared to the efficiency reached in 
Super-Kamiokande for the same decay channel ^| . 

The activity for the proton decay is given by the ex- 
pression: 

A = sNptm/r (8) 

where e = 0.65 is the efficiency already explained; Np = 
1.45 • lO'^^ is the number of protons in the detector; tm is 
the measuring time and r is the lifetime of the proton. 

For the current proton lifetime limit for the channel 
considered (r = 2.3-10^'^ y) 10], about 40.7 proton decay 
events would be observed in LENA after a measuring 
time of ten years with about 1.1 background events. If 
no signal is seen in the detector within this ten years, the 
lower limit for the lifetime of the proton will be placed at 
T > 4 • 10'^'* y at 90% C.L. If one candidate is observed, 
the lower limit will be reduced to r > 3 ■ 10'^'* y at 90% 
C.L. and the probability of this event being background 
would be 32%. 



VII. CONCLUSIONS 

Within the simulations performed, an efficiency of 
~ 65% for the search for proton decay in the LENA de- 



tector has been determined. A lower limit for the proton 
lifetime of r > 4 • lO^^ y (at 90% C.L.) can be reached if 
no proton decay event is measured within ten years. 

The success in reaching this limit is based on the dis- 
tinct pulse shape of the proton decay mode p — > K'^V 
which is clearly observable in the LENA detector. Only 
a minor energy cut was used and hence the efficiency for 
observing the proton decay including free and bound pro- 
tons is high. Since the values predicted by the favoured 
theories for the proton decay in this channel are of the 
order of the value resulting from our simulation, it is obvi- 
ous that LENA measurements would have a deep impact 
on the proton decay research field. In order to further 
improve the results of the present simulation pulse shape 
analysis of higher precision will be performed. 
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